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IV 


ABSTRACT 


Using  visible  wavelength  radiance  data  obtained  from 
the  spaceborne  Sea-viewing  Wide  Field-of-view  Sensor 
(SeaWiFS)  during  the  Aerosol  Characterization  Experiment- 
Asia  (ACE-Asia) ,  an  analysis  of  aerosol  optical  depth  (AOD) 
was  completed  by  modification  to  the  NFS  AOD  Model.  Cloud- 
free  AOD  estimates  have  been  compiled  in  previous  years 
from  NOAA  geosynchronous-  and  polar-orbiting  satellite  data 
and  validated  using  surface  sunphotometers .  The  objective 
of  this  thesis  was  to  calibrate  the  linearized,  single¬ 
scatter  algorithm  using  estimated  bi-directional  surface 
reflectance  and  size-dependent  phase  function  parameters. 
The  intent  of  the  study  was  to  provide  enhanced  temporal 
AOD  coverage  in  the  littoral  and  open  ocean  environment 
with  the  addition  of  the  orbiting  SeaWiES  eight -channel 
radiometer  to  the  established  NOAA  constellation  of  five- 
channel  AVHRR-equipped  satellites.  Comparison  to  Aeronet 
ground  stations  provides  in-situ  ground  truth.  "Clean" 
ACE-Asia  sky  regions  have  a  mode  at  SeaWiES  AOD  around 
0.25,  while  "dirty"  dust  plumes  had  a  mode  at  AOD  near  2.0, 
tailing  beyond  4.0.  Initial  SeaWiFS  AODs  were  about  20% 
higher  than  AVHRR  in  clean  subregions  and  up  to  100%  higher 
in  dirty  subregions.  Refined  ozone  and  Rayleigh  scatter 
parameters  have  reduced  SeaWiFS  excess  AOD  by  6%  to  12%. 
Red  tide  surface  effects  and  multiple  scatter  atmospheric 
effects  were  present,  complicating  current  assumptions. 

The  work  has  operational  significance  in  providing 
more  timely  remote  sensing  data  to  military  operators  of 
electro-optical  identification  and  targeting  systems. 
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I. 


INTRODUCTION 


Satellite  imagery  from  the  National  Oceanographic  and 
Atmospheric  Administration  (NOAA)  -16  and  Orbview-2  polar- 
orbiters  are  used  to  estimate  aerosol  optical  depth.  In 
this  thesis,  an  experiment  over  east  Asia  combines  the 
untested  SeaWiFS  sensor  algorithm  with  ground  station 
observations  and  the  well -calibrated  AVHRR/3  algorithm.  A 
first  iteration  of  corrections  to  the  SeaWiFS  product  are 
completed,  and  a  discussion  of  the  changing  weather 
conditions  and  variables  involved  in  this  process  follow. 

Crude  aerosol  detection  from  space  began  with  the 
first  successful  Suomi  radiometer  images  from  Explorer  7, 
launched  on  October  13,  1959.  Similar  direct  observations 
of  agricultural  burning,  desert  dust  storms,  and  volcanic 
eruptions  are  made  to  this  day.  Slightly  more  advanced 
limb  scattering  measurements  for  diffuse  stratospheric 
particles  were  first  hand-made  by  Apollo-Soyuz  astronauts 
in  1975,  and  the  initial  tropospheric  estimates  were  made 
using  ERTS-1  (Landsat)  the  same  year. 

While  precision  has  improved  and  spatial  resolution 
has  increased,  temporal  overhead  frequency  of  aerosol 
observations  is  still  limited  by  overhead  coverage. 
Geostationary  satellites  provide  whole-disc  coverage  every 
half-hour  at  the  expense  of  resolution  from  higher  altitude 
orbits.  Sunsynchronous  polar-orbiters  are  cheaper  to 
launch  and  thus  carry  the  latest,  highest -quality 
instruments.  Their  low-earth  orbit  profile  does  limit  them 
to  only  view  about  half  the  earth  twice  daily  and  sometimes 
have  gaps  in  their  coverage. 
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A. 


CURRENT  DETECTION  OF  AEROSOL 


Directly  sampling  the  constituents  of  the  atmosphere, 
including  aerosols,  would  require  climbing  the  tops  of 
mountains,  launching  balloons,  or  flying  instrumented 
airplanes  to  detect  aerosol  particles.  We  use  remote 
sensing  to  gain  timely,  global,  and  relatively  inexpensive 
observational  coverage. 

NOAA' s  newest  polar-orbiter ,  N-16,  (Figure  1)  went 
into  orbit  in  September,  2000  carrying  the  Advanced  Very 
High  Resolution  Radiometer  (AVHRR/3) .  It  is  capable  of 
1.09  km  resolution  in  five  channels,  two  broad  wavelength 
bands  of  which  are  partly  in  the  visible  light  spectra. 
She  orbits  at  870  km  above  earth  in  a  98.7°  retrograde 
orbit,  passing  the  equator  northward  (ascending  node)  at 
1:40  PM.  Her  precession  of  0.986°  per  day  keeps  the  sun  at 
a  constant  angle  in  the  west  on  her  afternoon  passes.  The 
AVHRR  views  ±  55.3°  from  nadir  for  a  swath  width  of  2,400 
km.  About  a  third  of  the  swath  is  affected  by  sunglint  to 
varying  degrees  depending  on  latitude  and  wind  speed. 
Storms  causing  large  wave  heights  increase  the  sunglint 
extent  due  to  reflection  off  individual  crest  facets. 


Figure  1.  NOAA-16  and  AVHRR/3.  (From:  Brown) 
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Orbital  Imaging  Corporation's  OrbView-2  (Figure  2)  was 
launched  in  August,  1997  carrying  the  Sea-viewing  Wide 
Field-of-view  Sensor  (SeaWiFS) .  Its  eight  narrow  visible- 
spectrum  wave  bands  provide  1.13  km  resolution  from  a  noon 
ascending  node  705  km  above  the  earth.  It  sees  a  swath 
2,801  km  wide  by  sweeping  ±  58.3°  over  a  point  ±20°  off 
nadir,  thus  avoiding  sunglint. 


Figure  2.  OrbView-2  and  SeaWiFS.  (From:  Hooker) 

These  space  sensors  must  be  calibrated  to  ground  truth 
to  assure  accuracy.  Local  measurements  are  made  under  a 
variety  of  conditions  -  marine,  desert,  forest,  or 
mountainous  terrain  -  depending  on  the  remote  sensing 
objective.  For  ACE-Asia,  the  AERONET  program  deployed  a 
network  of  sunphotometers  shown  in  Eigure  3  from  the  Gobi 
Desert  at  the  foothills  of  the  Himalayan  Mountains  eastward 
to  islands  in  the  western  Pacific  Ocean.  More  details  of 
these  sensors  will  be  discussed  in  Chapter  III. 
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Figure  3.  ACE-Asia  Network  and  Cimel  318-1. 


(From:  Bates  et.al.) 

B.  IMPORTANCE  OF  AEROSOL  DETECTION 

Aerosols,  in  addition  to  causing  health  and  pollution 
concerns,  act  indirectly  to  counter  the  "greenhouse" 
warming  effect  by  increasing  reflected  solar  energy  back  to 
space.  Though  tiny,  individual  particles  scatter  light 
according  to  their  size,  shape,  and  chemical  composition. 
From  high  overhead,  the  satellite  measures  solar  radiation 
scattered  out  of  the  atmosphere.  The  ocean  appears  dark, 
while  aerosols  are  a  spatially  and  temporally  variable 
source  of  scattered  radiance.  Aerosol  particles  are  kicked 
up  by  high  winds  or  lofted  from  anthropogenic  and  natural 
sources,  carried  across  thousands  of  miles  by  atmospheric 
currents,  and  ultimately  fall  out  of  the  air  column  in 
precipitation  and  dry  deposition  processes.  While  thick 
smoke  and  pollution  appears  dark  in  color  imagery,  sea 
salts  and  haze  brighten  images  by  increasing  backscatter 
towards  the  sensor  much  like  a  thin  cirrus  layer. 

Military  and  private-sector  interests  alike  are  more 
reliant  than  ever  on  remote  sensing  and  optical  systems  for 
weaponry,  communications,  safety,  and  logistics.  A  precise 
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knowledge  of  the  aerosol  conditions  over  foreign  shores  and 
above  our  own  heads  is  critical  for  saving  lives  and 
increasing  efficiency.  The  results  of  this  experiment 
might  have  lasting  effects  on  littoral  optical  forecasts 
for  the  warfighter  and  search-and-rescue  personnel,  or  for 
treaty  enforcement  between  governments. 
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II.  AEROSOL  OPTICAL  DEPTH  RETRIEVAL 


A.  RADIATIVE  TRANSFER  THEORY 

Energy  reaching  a  satellite's  radiometer  comes  from 
solar  reflection  or  terrestrial  emission.  Sunlight  enters 
the  earth' s  atmosphere  and  scatters  back  to  space  or 
reflects  off  the  planet's  surface.  At  longer  wavelengths, 
infrared  and  microwave  emission  from  the  ground  and 
atmospheric  constituents  are  important,  but  emission  is  not 
important  in  the  visible  region.  The  reflected  intensity 
changes  due  to  scatter  into  or  out  of  the  beam  as  well  as 
loss  along  the  path  due  to  absorption.  Little  sunlight  is 
reflected  from  the  ocean  surface  except  in  a  narrow 
sunglint  region  defined  by  the  geometry  required  for 
specular  (mirror-like)  reflection.  Like  the  greater 
reflection  from  clouds  and  land  surfaces,  the  small  diffuse 
portion  of  ocean  surface  reflectance  increases  the  expected 
brightness  level  received  at  the  sensor  beyond  what  would 
be  expected  from  the  gaseous  atmosphere  alone.  A  further 
increase  can  come  from  diffuse  multiple-scattering  from 
just  below  the  ocean  surface,  when  light  interacts  with 
suspended  particles  in  the  water  column.  As  will  be 
discussed  further,  this  complicates  littoral  minehunting 
operations  when  turbid  outflow  or  plankton  blooms  are 
present.  These  types  of  events  are  common  along  the  Asian 
coast,  giving  the  Yellow  Sea  its  name. 

Figure  4  indicates  with  colored  bands  the  eight 
SeaWlFS  channels  plotted  over  portions  of  the  atmosphere 
which  allow  high  energy  transmittance.  These  channels  in 
the  visible  portion  of  the  spectrum  are:  0.402-0.422, 


7 


0.433-0.453,  0.480-0.500,  0.500-0.520,  0.545-0.565,  and 
0.660-0.680  |lm.  In  the  near-infrared,  SeaWiFS  measures 
wider  bandwidths  at  0.745-0.785  |lm  and  0.845-0.885  |Jm. 


Figure  4.  SeaWiFS  channels  over  atmospheric  windows. 

Of  particular  interest  are  SeaWiFS  channels  six  (0.670 
|fm)  and  eight  (0.865  |im)  ,  shown  in  Figure  5,  which 
correspond  roughly  to  AVHRR  channels  one  (0.580-0.680  |lm,  a 
lower  frequency  and  five  times  wider  bandwidth)  and  channel 
two  (0 . 725-1 . 100|lm,  nearly  equal  center  wavelength  but 
nearly  100  times  wider  bandwidth) . 


Figure  5.  AVHRR  channels  over  solar  irradiance. 


8 


Band  differences  cause  sensors  to  receive  different 


amounts  of  reflected  solar  radiation,  as  shown  by  the 
representative  areas  under  a  theoretical  Plank  curve  in 
Figure  6. 


Figure  6.  Radiometer  radiance  reception  limited  by 

channel  width. 

Particle  size  plays  an  important  factor  in  light 
interaction  as  well.  Light  at  increasing  wavelengths 
interacts  with  particles  of  increasing  size  in  one  of  three 

regimes.  A  size  parameter  is  defined  as  particle  diameter 

divided  by  wavelength  (%  =  2Kr  I X  )  .  Scatter  from  cloud 

droplets  and  rain  falls  into  the  largest  visible  wavelength 
size  parameter  category,  geometric  optical  scatter.  This 
effect  explains  the  common  occurrence  of  rainbows  and  sky 
halos.  Atmospheric  molecules  and  tiny  Aitken  particles 
cause  the  smallest  size  parameter  category,  Rayleigh 

scatter.  In  between  lies  the  Mie  scatter  regime  in  which 
dust  and  smoke  particle  size  is  roughly  equal  to  the 
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wavelength  of  the  visible  light.  This  also  explains  why 
electromagnetic  radiation  at  longer  wavelengths,  such  as 
the  upper  infrared,  microwave,  and  radar  regions,  have 
little  interaction  with  aerosols. 

This  thesis  assumes  sunglint-free,  cloud-free,  clear 
water  conditions  and  spherical  aerosol  particles.  Constant 
surface  reflection  from  sea  foam  caused  by  high  winds  and 
phytoplankton  blooms  is  ignored.  Correction  factors  for 
Rayleigh  scatter  and  absorption  due  to  ozone  and  water 
vapor  are  applied. 


B.  AN  INVERSE  SOLUTION  TO  THE  RADIATIVE  TRANSFER  PROBLEM 


Aerosol  optical  depth  can  be  retrieved  through  an 
inverse  process  if  the  total  measured  radiance  (Lt)  is 
compared  to  the  known  solar  constant  (Eq)  . 

1 .  Aerosol  Radiance 

Total  extinction  is  the  sum  of  absorption  and 
scattering  (0^  =  a^+0^)  depending  on  wavelength.  Radiance 

decreases  proportional  to  the  original  intensity  along  the 
path  according  to  Beer's  law: 

UX)  = 

Optical  depth,  then,  is  the  integrated  extinction 
coefficient  for  all  wavelengths  along  the  path: 
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(2) 
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The  summation  of  radiation  received  by  the  satellite 
is  shown  in  Figure  7  as  the  aerosol  backscatter  of  interest 
(not  shown) ,  the  Rayleigh  molecular  scatter  (not  that  of 
light  from  the  surface  scattering  into  the  sensor's  field 
of  view) ,  and  the  percentage  of  surface  reflected  light  and 
sunglint  that  the  atmospheric  transmissivity  allows  to  pass 
through . 


Radiathe  Transfer  Theors’ 

% 


Figure  7 . 


A  summation  of  the  NFS  model  treatment  of  the 
radiative  transfer  theory.  (After:  Durkee) 
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Liou  (1980)  gave  the  general  form  of  total  radiance 


through  the  atmospheric  plane  as: 

f 


dL^(k,£l)  _  j  .r.  03^ 


jL,(^,n)/)(\|/gdn  +  E^%p(y/ ^)e  ^ 


4ji 


(3: 


where  ?o  is  the  single-scatter  albedo  (ratio  of  scatter 
extinction  to  total  extinction)  and  ignoring  multiple 
reflections  in  clouds  or  undersea.  P  (¥s)  is  the  scatter 
phase  function  for  the  single  scattering  angle  (\|/s)  between 
the  incoming  solar  irradiance  and  the  satellite. 


Solving  the  first  order,  linear,  ordinary  differential 
equation  for  aerosol  sources  of  airborne  scatter  only: 


L  ~  L 

c>  ,  ^  sun 

471 


{l-e 


-5(X,z)(1/|x„+1/H)^ 


(1/|I,  +1/|I) 


(4) 

as  notated  by  Durkee  et  al .  (1991) .  In  this  single-scatter 
limited  condition,  the  total  radiance  received  includes  in 
the  first  term,  surface  reflected  energy  not  scattered  out 
of  the  path.  The  second  term  includes  the  percentage  of 
light  that  has  scattered  back  from  the  atmosphere  and 
happens  to  go  in  the  direction  of  the  receiver  without 
transmitting  to  the  surface  or  being  absorbed. 
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Simplified,  if  the  aerosol  optical  depth  is  assumed  to  be 
much  less  than  one,  the  aerosol  radiance  becomes: 

^aerosol  .  ^ sun  (  S  ■) 

471  ILl  ^  ^ 

from  Kidder  &  Vender  Harr  (1995)  .  La  is  proportional  to  ?o, 
p  (Vj/s)  ,  and  da,  and  the  increased  brightness  measured  from 
space  is  accounted  for  by  the  optical  depth  of  the  aerosol 
modified  by  the  directionality  of  the  scatter  and  the 
particle  absorption. 

2 .  Backscatter  Parameterization 

Backscatter  towards  the  sensor  is  one  of  the  keys  to 
the  inverse  solution.  The  NFS  model  currently  assumes  only 
simple  "single-scattering"  of  light  off  of  spherical 
aerosol  and  back  towards  the  satellite  sensor.  The  amount 
of  energy  scattered  and  absorbed  off  an  aerosol  surface  is 
determined  by  the  index  of  refraction  (m  =  n  +  in').  In  the 

visible  and  near-infrared  portions  of  the  spectrum, 
absorption  is  negligible.  Scattering  includes  not  only 
photons  which  strike  a  particle,  but  those  which  refract 
around  while  passing  near  it.  Double  integrals  in  the 
solution  to  Eg.  (3),  even  neglecting  irregular  shapes,  are 
difficult  to  estimate  since  knowledge  of  radiance  into  the 
scattering  volume  from  all  directions  is  required. 
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One  way  to  get  around  this  problem  is  the  use  of  look¬ 
up  tables  for  phase  function  values  determined  by  a  variety 
of  aerosol  types  and  size  distributions.  Since  the 
scattering  is  determined  by  the  product  of  the  Mie 
scattering  efficiency  which  peaks  near  particle  radii  equal 
to  the  wavelength,  the  more  steeply  the  number  of  particles 
falls  off  as  size  increases,  the  greater  the  change  in 
scattering  as  a  function  of  wavelength.  Figure  8  shows  the 
process  of  determining  an  optical  depth  using  the  NFS 
model . 

NPS  Aerosol  Model 


orbit  parameters; 


Figure  8 .  Look-up  table  method  for  approximating  phase 

function.  (After:  Durkee) 
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The  ratio  of  AVHRR  channels  one  and  two  (or  SeaWiFS 
channels  six  and  eight)  determines  the  size  distribution 
and  phase  function  for  the  scatter  angle  determined  by  the 
satellite-solar  geometry.  This  estimated  phase  function 
allows  calculation  of  an  aerosol  optical  depth: 


1  1 

- ILl 

471 


(6) 


If  a  higher  size  distribution  number  is  selected,  the 
result  is  that  the  model  assumes  a  larger  percentage  of 
small  particles  are  present  in  the  atmosphere.  These 
smaller  particles  have  the  size  parameter  (%)  closer  to 
unity.  The  increased  phase  function  assumed  by  the  model 
on  the  right  side  of  the  bottom-left  plot  would  increase 
for  values  of  scattering  angle  over  140°.  Located  in  the 
denominator  of  Equation  6,  the  resulting  aerosol  optical 
depth  calculation  would  be  decreased.  Conversely,  at 
scattering  angles  less  than  140°,  an  error  in  increased  size 
distribution  would  decrease  phase  function  and  increase  the 
output  AOD . 

A  "two-stream  Turner"  routine  accounts  for  ozone  and 
Rayleigh  scatter  contributions  to  radiance  measured  at  the 
receiver.  A  constant  Rayleigh  phase  function  formula  of 
^^(l  +  cos  2\|/s)  is  used  in  this  calculation.  For  a  more 
detailed  treatment  of  the  mult ispectral  parameterization 
routine  used  in  the  NFS  model,  see  Brown,  98,  pp .  24-27. 
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Ill .  EXPERIMENTAL  DATA 


A.  PURPOSE  OF  ACE -ASIA 

Global  aerosol  chemistry  and  transport  measurements 
have  been  undertaken  by  the  International  Global 
Atmospheric  Chemistry  Program.  Previous  investigations 
have  focused  on  ocean  regions  around  Tasmania  (ACE-I),  the 
Canary  Islands  (ACE-II),  the  east  coast  of  the  United 
States  (TAREOX) ,  and  the  Indian  Ocean  (INDOEX).  ACE -ASIA 
was  held  from  March  to  May  2001  in  order  to  quantify 
aerosol  properties  and  their  radiative  effects  over  East 
Asia  and  the  western  North  Pacific  Ocean. 

The  winter  Himalayan  high  drives  northwest  gales 
across  the  Gobi  Desert.  Tons  of  fine  western  desert 
mineral  sand  is  joined  by  coal  soot  and  SO2  from  the 
agricultural  heart  of  China.  finally,  the  current  picks  up 
NO2  emissions  in  the  form  of  pollution  from  industrial 
factories  and  the  auto  traffic  of  crowded  eastern  Asian 
cities.  Weather  patterns  carry  these  dirty  air  masses  out 
over  the  Yellow  Sea,  the  Sea  of  Japan,  and  several  days 
later  can  even  cross  the  Pacific  to  the  west  coast  of  the 
United  States,  as  seen  in  figure  9. 
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Figure  9.  Dust  storm  crossing  the  Pacific  Ocean. 

Sea  of  Japan,  10  APR  01;  Southern  California,  16  APR  01. 

(From:  Hooker) 

An  unusual  AOD  observation  of  the  ACE-ASIA  experiment, 
besides  the  high  AOD  values,  was  the  four-five  day  period 
of  dust  plumes  following  a  frontal  passage  of  the  eastern 
Aeronet  sites.  Also,  the  western  Aeronet  sites  experienced 
diurnal  shifts  from  morning  coal  soot  and  industrial 
sulfates  to  afternoon  and  evening  dust  concentrations. 
Since  this  thesis  includes  only  local  afternoon  overpasses, 
a  mixture  of  the  two  particle  types  can  be  expected. 
Aerosol  visible  single-scatter  albedo  differences  between 
African  Saharan  dust  with  sea-salt  and  nitrate  mixes  ((Oq  = 
0.8  to  0.9  during  ACE  II  by  Collins,  et .  al .  )  and  Asiatic 
Gobi  Desert  dust  with  industrial  sulfate  and  coal  soot 
mixes  (from  0.3  to  0.95  by  Xu  and  Bergin,  2001)  have  been 
reported.  The  effect  of  particle  albedo  in  a  visible  light 
Mie  scatter  regime  needs  to  be  further  investigated. 
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B .  INSTRUMENTS 

ACE-Asia  joined  scientific  minds  and  resources  from 
over  ten  countries.  Terrestrial  data  collection  was 

accomplished  from  land,  sea,  and  air.  In  addition  to  the 
ground  station  radiometers  and  LIDARs,  the  R/V  Ronald  H. 
Brown  and  R/V  Mural  were  involved,  as  well  as  aircraft 
including  the  NCAR  C-130,  NFS  CIRPAS  Twin  Otter,  and 
Australian  King  Air. 

1.  NOAA  AVHRR/3 

Channels  one  (0.580  -  0.680  |lm)  and  two  (0.725  -  1.100 
|j,m)  of  the  AVHRR/3  are  used  to  estimate  aerosol  optical 
depth.  A  split-window  cloud  mask  routine  uses  channels 
one,  two,  four  (10.3  -  11.3  |lm)  ,  and  five  (11.5  -  12.5  |lm)  . 
The  spectral  response  curves  for  channels  one  and  two  are 
shown  in  Figure  10. 


500  550  «00  650  7C0  750  SOO 

vviwelErt^li  (N  anometef i) 


6E0  740  SOO  560  P20  980  3040 

(Nancmctert) 


Figure  10.  AVHRR  spectral  response  (chs.  1  &  2)  . 

(From:  Kidwell) 
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Experimental  analysis  yields  an  ozone 
optical  depth  0.027  in  channel  one  and  0.0021 
two.  The  effects  of  Rayleigh  scatter  optical 
assumed  to  be  0.057  in  channel  one  and  0.019 
two . 


absorption 
in  channel 
depths  are 
in  channel 


2 .  NASA  SeaWiFS 

A  successor  to  the  Coastal  Zone  Color  Scanner  on 
NIMBUS-7,  this  multi-band  color  monitor  can  track 
meteorological  events  as  well  as  detect  subtle  differences 
in  biological  activity  in  the  world's  oceans.  Actually 
built  and  owned  by  Orbital  Sciences  Corporation,  NASA  buys 
and  owns  the  right  to  all  research  data  collected. 

In  order  to  extend  coverage  of  bright  areas  without 
saturating  the  scanner,  SeaWiFS  uses  discontinuous  gain 
above  80%  peak  input.  This  bi-linear  response  allows,  for 
example,  channel  one  to  detect  upper  tropospheric 
cloudtops,  ice  fields,  or  desert  sands  at  up  to  60.1  mW  cm^^ 
pm^^  sr^^  vice  the  original  limit  of  about  13.6  mW  cm^^  pm^^ 
sr^^.  The  spectral  response  curves  for  SeaWiFS  channels  six 
and  eight  are  shown  in  Figure  11  for  comparison  with  AVHRR. 
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Figure  1 1 . 


SeaWiFS  bands  six  and  eight. 
(Author) 


The  ability  of  SeaWiFS  to  avoid  sunglint  by  automated 
tilting  is  an  advantage  over  AVHRR,  but  it  can  cause 
unusual  streaks  in  SeaWiFS  imagery  as  the  tilt  angle  is 
changed . 

3.  Aeronet  Cimel  CE  318-1 

The  Aerosol  Robotic  Network  (Aeronet)  radiometer 
instrument  detects  transmitted  solar  intensity  in  five 
visible  and  three  near-infrared  channels  (0.340,  0.378, 

0.440,  0.499,  0.613,  0.870,  0.940,  and  1.020  pm) 

(ARM/NASA/DOE  website,  after  Holben,  1997).  AOD  estimates 
are  possible  with  an  accuracy  of  -i—  0.02.  Aerosol  modal 
radius  distribution  information  from  0.1  to  3  pm  is 
estimated  using  a  simple  two-step  radiative  transfer 
algorithm  which  had  not  been  fully  validated  as  of  June, 
2001 

Each  unit  uses  a  four-quadrant  feedback  system  to 
center  on  the  sun  every  fifteen  minutes  during  daylight 
hours  and  sample  incoming  solar  radiance  using  two 
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polarized  collimator  filters  in  front  of  the  sensor.  The 
original  Aeronet  sites  in  the  continental  United  States  are 
funded  by  the  Department  of  Energy' s  Atmospheric  Radiation 
Measurement  program.  These  sensors  have  the  ability  to 
transmit  their  readings  automatically  over  an  uplink  to  a 
NOAA  GOES  relay.  Nine  remote,  manned  sites  were  combined 
with  local  data  collection  teams  during  ACE-Asia. 

Data  output  screening  is  performed  using  a  number  of 
quality  checks.  Eirst,  samples  are  always  taken  in 
"triplets"  separated  by  thirty  seconds.  Any  sudden  changes 
in  AOD  indicates  cloud  contamination.  AODs  less  than  0.01 
and  very  low  sun  angle  measurements  near  twilight  are 
thrown  out  to  avoid  the  associated  variability  which  would 
bias  daily  averages.  Results  greater  than  three  standard 
deviations  from  the  diurnal  mean  are  discarded.  However 
meticulous,  these  standards  can  allow  thin,  stable,  uniform 
clouds  to  pass  and  misrepresent  thick,  quickly  arising  dust 
columns  as  cloud  streaks. 
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IV.  RESULTS 


A.  IMAGERY  RETRIEVAL  AND  PROCESSING 

Data  collected  during  ACE-Asia  and  used  in  this  study 
includes  reference  ground  data  sets,  raw  satellite  radiance 
imagery,  and  supporting  meteorological  and  oceanographic 
fields.  The  "ground  truth"  state  of  the  atmosphere  and 
underlying  ocean  are  needed  for  accurate  validation  of 
aerosol  optical  depth  retrievals.  Ocean  color  variability 
affects  aerosol  measurements  at  the  same  time  that  aerosol 
changes  affect  water  color  data  collection. 

Preliminary  data  processing  was  accomplished  with  the 
Terascan  software  suite.  For  AVHRR  and  SeaWiFS  images,  the 
initial  step  was  to  build  information  on  the  swath  from  sun 
and  satellite  angles  at  each  location  being  recorded. 
Precise  spatial  knowledge  is  required  for  image 
intercomparison.  Radiance  to  reflectance  conversion  is 
applied  based  on  wavelength  dependence  on  Eq,  day  of  year 
corrections  for  annual  variation  in  incoming  solar 
intensity,  and  satellite  view  geometry  is  used  to 
georeference  the  image. 

B.  NPS  MODEL  CODE  PARAMETERIZATION 

Throughout  the  study,  data  quality  was  favored  over 
quantity,  and  sections  or  entire  passes  were  removed  to 
minimize  erroneous  analysis.  Sunglint  removal  was 
completed  to  remove  anisotropic  specular  reflection  into 
the  receiver.  Cloud  screening  was  completed  by  a 
combination  of  automated  and  manual  techniques.  Due  to  the 
lack  of  infrared  bandwidth  receivers  on  SeaWiFS,  a  split- 
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channel  method  was  partially  ineffective.  Previous  AVHRR 
routines  have  compared  cloudtop  to  surface  temperatures,  at 
the  risk  of  confusing  stratus  and  fog.  Another  method  has 
been  spatial  coherence  analysis  to  find  cloudtop  texture, 
though  this  might  mistake  thick  aerosol  or  turbid  sea 
surface  filaments. 

A  visual  analysis  showed  that  patches  of  low  closed¬ 
cell  cumulus  in  the  open  ocean  to  the  south  and  east  of 
Japan  often  passed  the  automatic  filter.  Also,  high 
turbidity  from  the  Yangtze  River  carried  into  the  western 
Yellow  Sea  renders  composite  AODs  meaningless  due  to 
excessively  high  surface  reflectance.  Manual  thresholds 
were  required  to  remove  bad  pixels,  even  at  the  expense  of 
some  good  data.  In  the  end,  average  "clean"  skies  had  AODs 
below  about  0.25,  "medium"  periods  had  AODs  up  to  0.50,  and 
"dirty"  events  had  AODs  as  high  as  4.0  and  above.  These 
extremes  were  corroborated  by  the  AERONET  ground  stations. 
In  comparison,  on  a  hazy  day  over  a  large  city  on  the  east 
coast  of  the  United  States  the  AOD  generally  doesn't  rise 
above  about  0.5.  All  products  included  in  this  thesis  have 
been  limited  to  an  AOD  of  2.0  in  order  to  show  finer  detail 
in  weaker  aerosol  regions. 

C .  CRITICAL  EVENTS 

More  than  120  SeaWiES  and  NOAA-16  overpasses  were 
downloaded  during  45  days  of  ACE -ASIA  operations.  Imagery 
had  to  be  selected  to  best  represent  the  range  of  aerosol 
observed.  This  study  required  cloud-free  passes  which 
coincided  in  swath  coverage  and  proximity  to  ground 
stations.  Many  swaths  did  not  overlap  with  each  other,  and 


often  ground  sites  could  not  report  readings  due  to  local 
cloudiness  that  did  not  extend  far  across  the  coastline. 
Unfortunately,  the  same  remote  sites  that  were  favorably 
far  from  industry  and  population  centers  frequently  were 
affected  by  orographic  clouds  due  to  their  hilltop 
locations.  Twenty-seven  matches  were  found  between  the 
sensors  and  these  were  narrowed  to  five  choices  that  had 
the  added  benefit  of  closely  matching  overpass  times, 
thereby  avoiding  advection  of  small  aerosol  features 
between  images. 

Two  of  these  aerosol  events  have  been  selected  to 
indicate  the  NFS  AOD  model  performance  under  significantly 
different  conditions.  A  third  sample  shows  an  extreme  dust 
plume  event.  Prevailing  conditions  during  the  spring 
transition  from  winter  to  summer  monsoon  behavior  are  seen 
in  Figures  12  &  13.  During  the  45-day  experiment  period, 
cloud  patterns  show  frequent  lows  developing  to  the 
northeast  of  the  Tibetan  highlands,  proceeding  to  sea  over 
Korea  and  Japan.  To  the  south  and  east  of  Japan  the 
conditions  are  opposite:  easterlies  and  surface  highs  form 
generally  clear  skies  except  for  cumulus  cells  drifting 
from  the  vicinity  of  Okinawa  towards  Taiwan.  Dust 
outbreaks  generally  follow  these  northern  lows,  trailing 
their  associated  cold  fronts  toward  the  southwest. 
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Figure  12 . 


"Typical"  cloud  patterns:  GMS  April  5,  2001 


Topographic  &  Mesoscale  Effects 


Figure  13.  Prevailing  spring  conditions  combined  with 

topographic  effect,  left;  mesoscale  wind  pattern,  right. 
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1. 


Similar  AOD,  Different  Look  Angle:  Day  115. 


The  first  day  selected  shows  similar  AOD  values  and 
the  effect  of  increased  column  distance  on  remote  sensing. 
The  location  of  interest  for  this  case  is  37°  37'  North 

latitude,  136°  36'  East  longitude.  On  Figures  14  and  15, 
the  point  chosen  for  comparison  is  signified  by  the  red 
star.  The  scattering  angle  is  128°  on  the  SeaWiFS  pass, 
while  NOAA-16  has  a  scattering  angle  of  157° 


Figure  14.  NOAA-16  pass  20011150433Z;  April  25,  2001. 
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Figure  1 5 . 


SeaWiFS  pass  20011150433Z; 


April  25, 


2001. 
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a.  Weather  at  Nearest  Aeronet  Station 

Surface  NOGAPS  fields  for  case  one  are  shown  in 

Figure  16. 


25  April:  NOGAPS  OOZ  Pro^ 


Figure  16.  Meteorological  conditions  for  April  15,  2001. 

Unfortunately,  although  weather  conditions  had 
been  calm  for  several  previous  days,  only  Noto,  Japan  had 
AERONET  data  nearby.  This  distance  from  Noto  to  the  swath 
overlap  was  at  least  375  kilometers,  upwind,  and  so 
comparison  between  the  two  are  marginally  useful.  Noto 
values  peak  at  0.28,  about  25%  lower  than  either  satellite. 
Travel  at  the  rate  of  this  wind  speed  would  take  more  than 
72  hours.  This  is  obviously  too  long  for  steady-state  air 
advect ion . 
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b.  NPS  Model  Output 

Aerosol  optical  depth  output  for  SeaWlFS  and 
AVHRR  are  in  the  top  row  of  Figure  17.  On  the  bottom  left 
is  the  difference  between  the  two,  and  on  the  bottom  right 
is  a  corresponding  GMS  image  showing  regions  masked  by 
cloud . 


25  April;  WiFS:  0433Z,  N16:  0433Z 


Figure  17.  SeaWiFS  AOD,  NOAA-16  AOD,  AOD  difference,  and 
GMS  showing  cloud  cover  for  April  25,  2001. 

This  set  of  AVHRR  and  SeaWiFS  passes  are  unique 
for  the  study  in  that  they  have  exactly  the  same  overpass 
time.  Generally  low  AOD  values  in  light  blue  on  the  top 
two  images  are  observed  over  the  eastern  Sea  of  Japan. 
Initial  SeaWiFS  parameterization  coefficients  produced  an 
AOD  value  (top  left  plot)  of  0.50;  improvements  to  ozone 
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and  Rayleigh  scatter  parameterization  brought  this  down  to 
0.46,  a  decrease  of  8.0%.  The  AVHRR  optical  depth  value 
(top  right  plot)  of  0.44  was  4.5%  lower  still.  These 
tendencies  were  common  throughout  the  entire  experiment, 
with  SeaWiFS  output  greater  by  up  to  20%  higher  over  clean 
skies  and  open  water,  becoming  worse  with  higher  aerosol 
optical  depths.  The  difference  in  these  AOD  outputs  is 
shown  on  the  bottom  left,  showing  low  excess  error  from 
SeaWiFS  in  low  dust  environments,  as  expected. 

The  bottom  left  plot  from  Figure  18  is  plotted 
again  in  Figure  18,  this  time  in  relative  percentage 
difference  above  AVHRR  instead  of  in  absolute  AOD  amount. 
Anomalous  high  values  south  of  Pyongyang  (#1)  are  caused  by 
turbidity,  Vladivostok  (#2)  caused  by  dust  and  cloud  edge. 


Figure  18.  Percent  difference  between  SeaWiFS  and  AVHRR 

AOD  for  April  25,  2001. 
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Refined  parameters  for  Rayleigh  scatter  and  ozone 
extinction  reduced  the  excess  SeaWiFS  AOD  values  by  6%  to 
12%.  Figure  19  shows  the  latest  iteration  of 
parameterization  with  AOD  coefficients  applied  for  this 
region  of  interest.  The  divisions  are  due  to  smoothing  on 
the  part  of  the  algorithm:  large  gradients  between  adjacent 
pixels  are  mistaken  by  the  computer  as  cloud  edge,  thus 
they  are  deleted  from  the  image. 


Figure  19. 


Latest  iteration  of  NFS  AOD  algorithm  for 
April  25, .2001 
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The  SeaWiFS  channel  six  value  for  this  case  is 
5.6  W/m^  sr  m,  channel  eight  is  2.49  W/m^  sr  m.  The  ratio 
of  channels  six  and  eight,  known  as  S12,  is  1.59.  The  phase 
function  therefore  is  0.15  using  a  backscatter  angle  of  128° 
and  size  distribution  model  number  5.64.  With  total 
radiance  of  60.18  W/m^  sr  m,  Rayleigh  radiance  removed  as 
19.45  W/m^  sr  m,  and  the  remaining  aerosol  radiance  is  40.72 
W/m^  sr  m.  Calculating  Equation  6  gives  an  AOD  of  0.49  for 
SeaWiFS.  NOAA-16  channel  one  radiance  value  is  4.69  W/m^  sr 
m,  channel  two  is  2.66  W/m^  sr  m  ,  for  an  S12  ratio  of  1.34. 
For  a  total  radiance  21.59,  Rayleigh  radiance  9.29,  the 
aerosol  radiance  is  10.51  for  phase  function  look-up  of 
0.26.  The  output  AOD  is  0.44,  having  selected  model  number 
1  with  backscatter  angle  157°. 

A  "light"  dust  day,  AOD  error  was  small  in  this 
case.  The  difference  in  retrieved  AOD  is  small,  even 
though  the  aerosol  radiance  values  are  quite  different. 
This  is  caused  by  entering  the  S12  lookup  tables  at  widely 
separated  look  angles  and  returning  nearly  identical  phase 
functions  from  two  different  modal  radius  curves  which 
balance  each  other  in  the  end. 

2.  Similar  Look  Angle,  Different  AOD:  Day  102. 


This  second  event  shows  closely-spaced  overhead  passes 
near  the  same  location  to  isolate  the  potential  causes  for 
differences  in  AOD  output  between  the  two  sensors.  The 
location  of  interest  for  this  case  is  33°  15'  North 
latitude,  126°  54'  East  longitude.  Cheju  Island,  is  seen  in 
Eigures  20  &  21.  Located  south  of  Punsan,  South  Korea  and 
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west  of  Sasebo  and  Nagasaki,  Japan,  the  site  was  ideal  in 
its  central  location  and  as  a  ground  sunphotometer  site. 
However,  it  had  a  unique  impediment  to  AOD  retrieval,  as 
well.  This  was  the  first  time  during  an  ACE  experiment 
that  scientists  encountered  the  well-known  possibility  of 
surface  reflection  by  phytoplankton  blooms. 


Figure  20 


NOAA-16  pass  20011020508Z,  April  12,  2001 
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Red  tides  were  seen  to  have  been  flowing  from  the 
Chinese  coast  to  the  southeast  towards  the  Tsushima  Straits 
(Figure  22) .  A  yearlong  deep  flow  is  formed  by  the  Yellow 
Sea  Warm  Current  and  the  China  Coastal  Current.  Flowing 
counter  to  the  counter-clockwise  gyre  pattern  might  seem 
contrary  to  a  physical  oceanographer,  but  these  blooms  are 
quite  shallow  and  drift  at  the  will  of  the  prevailing  winds 
rather  than  ocean-basin  thermohaline  forcing. 
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Figure  22.  Red  tides  in  the  Yellow  Sea  from  the  R/V  Ron 
H.  Brown.  (Photo  courtesy  of  Piotr  Flatau) 

These  red  tides  raise  surface  reflectance,  especially 
at  red  visible  wavelengths.  Seawater  absorbs  the  near- 
infrared  wavelengths  so  well  that  it  can  generally  be 
thought  of  as  nonref lect ive  in  the  absence  of  other 
materials.  Phytoplankton  pigment  may  create  a  greater 
effect  in  SeaWiFS  channel  six  than  AVHRR  channel  one  due  to 
its  narrow  bandwidth  centered  at  0.670  |lm.  Complications 
arise  when  blooms,  oil  slicks,  foam  whipped  up  by  winds 
over  14  meters  per  second,  and  partially  dissolved  sediment 
are  present.  Regional  studies  in  this  area  are  hampered  by 
the  same  conditions  that  give  the  Yellow  Sea  it's  name-  up 
to  1.4  billion  metric  tons  per  year  of  sediment  transport 
from  the  Yangtze  River  alone. 
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Red  tides  often  bloom  in  24-48  hours  and  last  one  week 
or  longer,  often  dying  out  just  as  quickly  due  to  lack  of 
nutrients.  Figure  23  shows  an  interesting  three-day  series 
over  Cheju  Island  from  April  12,  2001  to  April  14,  2001. 


Figure  23.  Red  tides  and  coastal  turbidity  variation 

from  space;  April  12-14,  2001. 

Image  one  shows  the  largest  plankton  blooms  as  well  as 
aerosol  flowing  from  north  to  south  along  the  left  edge  of 
the  image.  Image  two  shows  less  plankton  and  an  obvious 
wind  shift  from  the  southwest.  Image  three  shows  even  less 
plankton  but  greater  turbidity  along  the  extreme  left  of 
the  image.  Decaying  plankton  change  color  (reflecting 
light  at  different  wavelengths)  from  green  or  red  to  yellow 
and  brown.  Sample  imagery  from  a  three-day  period  show 
changes  in  underlying  seawater  constituents  that  could 
carelessly  be  taken  for  aerosol. 
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a.  Weather  at  Nearest  Aeronet  Station 

The  winds  over  Cheju  on  April  12,  2001  were 
common  over  the  entire  experiment  period.  Seen  in  Figure 
24,  north-northwest  winds  vary  in  intensity  and  back  from 
the  west  with  the  passage  of  numerous  lows.  These  nearly 
always  passed  to  the  north  with  their  associated  cold 
fronts  lingering  a  day  or  two  behind.  The  Himalayan  high 
has  brought  clear  skies  on  this  day  and  good  conditions  for 
AOD  retrievals. 


Figure  24 . 


Meteorological  conditions  for  April  12,  2001. 
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b. 


NPS  Model  Output 


Figure  25  is  unique  for  the  study  in  that  they 
have  exactly  the  same  scattering  angles,  162°.  Higher  dust 
values  over  the  western  Sea  of  Japan  and  Yellow  Sea  are 
observed.  Initial  SeaWlFS  parameterization  coefficients 
produced  an  AOD  value  of  0.66  (upper  left);  modifications 
brought  this  down  to  0.58,  a  decrease  of  12.1%.  The  AVHRR 
optical  depth  value  (upper  right)  of  0.44  was  31.8%  lower 
still.  This  excessive  SeaWiFS  optical  depth  at  higher  AODs 
was  common. 


Figure  2 5 . 


SeaWiFS  AOD,  NOAA-16  AOD,  AOD  difference,  and 
Aeronet  data  for  April  12,  2001. 


Further  complicating  this  retrieval  were  the 
underlying  red  tide  filaments  seen  in  the  upper  left  of 
Figure  25,  and  more  clearly  as  blotches  in  the  center  of 
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Figure  26.  Compare  this  to  the  uniform  difference  values 
for  the  clean  skies  in  case  one. 


Figure  26. 


Percent  difference  between  SeaWiFS  and  AVHRR 
AOD  for  April  12,  2001. 


Differences  in  SeaWiFS  and  AVHRR  AOD  are  clearly 
seen  in  the  above  Figure  27  over  28,  after  the  improved 
Rayleigh  and  ozone  coefficients  have  been  applied. 
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Figure  2 7 . 


NOAA-16  AOD  product,  April  12,  2001. 


Figure  28.  Latest  iteration  of  NFS  AOD  algorithm  for 

April  12 , .2001 
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The  SeaWiFS  channel  six  value  for  this  case  is 
2.42  W/m^  sr  m,  channel  eight  is  0.95  W/m^  sr  m.  The  ratio 
of  channels  six  and  eight,  S12,  is  1.66.  The  phase  function 
thus  is  0.18  using  a  backscatter  angle  of  162°  and  size 
distribution  model  number  6.51.  With  total  radiance  of 
only  24.94  this  time,  Rayleigh  radiance  at  8.53,  and  the 
remaining  aerosol  radiance  16.41.  AOD  becomes  0.58  for 
SeaWiFS.  NOAA-16  channel  one  radiance  value  is  4.57  W/m^  sr 
m,  channel  two  is  2.72  W/m^  sr  m,  for  an  S12  ratio  of  1.27. 
For  a  total  radiance  21.08,  Rayleigh  radiance  8.38,  the 
aerosol  radiance  is  10.90  for  phase  function  look-up  of 
0.28.  The  output  AOD  is  0.44,  again  with  model  number  1, 
backscatter  angle  also  162°.  SeaWiFS  error  is  greater  as 
AOD  increases,  however  the  increased  surface  reflection 
contamination  probably  has  the  greatest  effect  on  AOD  error 
in  this  region. 

3.  Extremely  high  AOD  event:  Day  100. 

Compare  the  preceding  low  AOD  values  to  a  particularly 
dirty  event  on  Julian  day  100  shown  in  Figures  29  and  30. 
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Figure  30.  SeaWiFS  pass  20011000345Z;  April  10,  2001. 

The  location  of  comparison  in  this  case  (37°  32'  N, 
141°  18'  E,  is  off  the  eastern  coast  of  Honshu  northeast  of 
Tokyo  Bay.  The  majority  of  the  dust  plume  lags  behind  a 
frontal  passage  over  Japan,  with  a  narrow  leader  of  dust  to 
the  east.  Further  zooming  in  these  images  shows  cloud 
shadow  on  the  underlying  dust  plume,  fixing  its  altitude 
between  the  marine  boundary  layer  and  the  cloudbase.  A 
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glitch  in  the  image  downlink  can  be  seen  in  Figure  30 
through  North  Korea. 

a .  Weather  at  Nearest  Aeronet  Station 

Upper-level  westerlies  carried  the  aerosol  to  the 
east  across  Japan  above  the  NOGAPS  surface  wind  fields  in 
Figure  3 1 . 


Figure  3 1 . 


Meteorological  conditions  for  April  10,  2001. 


45 


b.  NPS  Model  Output 

Aerosol  optical  depths  in  Figure  32  have  been 
clipped  at  2.0,  but  actually  exceeded  4.5.  The  AOD  output 
was  lost  when  it  began  merging  with  clouds  associated  with 
a  low  forming  northeast  of  Vladivostok. 


Figure  32.  SeaWiFS  AOD  20011000345Z;  April  10,  2001. 

This  very  dirty  case  SeaWiFS  channel  six  value  is 
6.05  W/m^  sr  m,  channel  eight  is  2.91  W/m^  sr  m.  S12  is 
1.73.  The  phase  function  is  0.14  using  a  backscatter  angle 
of  130°  and  size  distribution  model  number  6.26.  With  total 
radiance  as  66.32,  Rayleigh  radiance  removed  as  10.88,  the 
remaining  aerosol  radiance  is  now  as  high  as  55.44. 
Equation  6  gives  a  high  AOD  of  1.47  for  SeaWiFS.  NOAA-16 
channel  one  radiance  value  is  6.33  W/m^  sr  m,  channel  two  is 
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4.33  W/m^  sr  m,  for  an  S12  ratio  of  1.27.  For  a  total 
radiance  28.13,  Rayleigh  radiance  6.67,  the  aerosol 
radiance  is  19.78  for  phase  function  look-up  of  0.23.  The 
AOD  output  is  only  0.799,  having  selected  low  model  number 
one  and  backscatter  angle  at  the  phase  function  curve 
inflection  point  near  142°.  SeaWiFS  error  appears  to 
increase  nonlinearly  as  AODS  continue  to  rise.  Multiple 
scatter  in  the  air  column  would  explain  this  AOD  retrieval 
output.  As  the  number  of  Mie  scatter  particles  in  a  given 
atmospheric  volume  increases,  the  chance  of  repeated 
scatter  towards  the  sensor  rises.  This  in  turn  increases 
the  total  radiance  observed,  leading  to  excessive  AOD 
output . 

D .  VALIDATION 

The  NFS  model  with  SeaWiFS  input  showed  some  skill  at 
measuring  aerosol  optical  depth  compared  to  AVHRR  and 
Aeronet  sunphotometer .  Iterations  of  the  corrected 

validation  process  will  continue  to  improve  the  Rayleigh 
and  ozone  parameterization.  Table  1  shows  the  latest 
values  used  in  this  thesis. 

Table  1.  Corrected  parameterization  for  NFS  SeaWiFS  AOD 

algorithm. 


N-16: 

Ch.  1 

Ch.  2 

,  P. 

().(K)5 

0 

Eo 

1628 

1030 

0.057 

0.019 

0.027 

0.002 

WiFS: 

Ch.  6 

Ch.  8 

P.. 

0.005 

0 

Eo 

1584 

10.30 

0.044 

0.019 

^03 

0.013 

0.002 
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Figure  33  shows  an  initial  scatter  plot  of  the  27 
matching  SeaWiFS  and  AVHRR  overpasses  prior  to  selecting 
the  three  critical  events.  It  was  assumed  that  SeaWiFS 
AODs  would  be  a  linear  factor  higher  than  AVHRR.  Simple 
parameterization  corrections  could  then  be  used  to  lower 
the  linear  best-fit  AOD  distribution  to  match  Aeronet  and 
AVHRR  AOD  values. 


Aeronet  Channel  Eight  (870  nm) 


Aeronet  AOD  Output 

Figure  33.  Initial  summary  of  SeaWiFS  and  NOAA-16  AOD 

compared  to  Aeronet  ground  Observation  at  870  nm,  with 
linear  best-fit  superimposed. 
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Figures  34  &  35  show  improved  parameterization  AOD 

(green  circles)  plotted  under  initial  SeaWiFS  AOD  values 
(blue  squares)  against  Aeronet  channel  six  and  channel 
eight  AOD.  Only  five  of  the  cases  with  Aeronet  data 
available  from  the  original  27  matches  are  plotted  for 
clarity.  Taiwan  (triangles)  was  consistently  plankton 
bloom-free  and  dust-free  and  had  the  lowest  error  of  any 
matching  comparison  and  negligible  change  due  to  improved 
parameterization . 


NPS  Model  vs.  Aeronet:  Ch  6 
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Figure  34 
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NPS  Model  vs.  Aeronet:  Ch  8 
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Figure  35.  SeaWiFS  AOD  comparison  with  Aeronet  ground 

observations  at  870  nm. 


SeaWiFS  error  does  not  appear  do  be  linear  with  increasing 
AOD  under  the  spherical  aerosol,  single-scatter  assumption. 

Table  2  lists  columns  of:  SeaWiFS  and  AVHRR  overpass 
time,  elapsed  time  separation,  nearest  Aeronet  Station, 
channel  six  and  eight  Aeronet  AOD  for  overpass  time, 
SeaWiFS  and  AVHRR  AODS,  latest  parameterization  SeaWiFS 
AOD,  and  relative  decrease  in  AOD  error  from  AVHRR. 
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Table  2  . 


Reduction  in  SeaWiFS  AOD  error  with  improved 
parameterization . 


WiFS 

N-16 

dT 

Aero 

AW6 

AN6 

AW8 

AN8 

W1 

m 

AOD 

DroD 

100  0345 
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D 
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0.421 
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1.39 

6.08% 

102  0333 
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95 

Cheju 
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0.262 

0.66 

0.44 

0.58 

12.12% 
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42 

Cheju 
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1.28 

0.76 

1.20 

6.25% 

108  0433 

108  0405 

-28 
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V. 


CONCUSIONS  &  RECOIMENDATIONS 


A  critical  analysis  of  the  modification  to  the  NPS 
aerosol  optical  depth  retrieval  model  for  SeaWlFS  input 
shows  promising  potential  with  a  few  shortcomings.  SeaWiFS 
AOD  values  were  higher  by  20-50%,  in  a  non-linear  manner 
at  high  AOD.  Future  construction  of  a  multi-dimensional  S12 
lookup  table  for  multiple  scatter  may  be  required  to 
correct  this  error. 

Size  distribution  model  numbers  were  larger  in  all 
cases,  as  well  as  S12  ratios.  Aerosol  radiance  for  case  one 
was  four  times  higher  for  SeaWiFS  and  50.6%  greater  for 
case  two.  These  factors,  along  with  backscatter  phase 
angle,  determine  the  phase  function  intrinsic  to  aerosol 
optical  depth  from  equation  six. 

Red  tides  combined  with  highly  reflective  turbid 
seawater  complicates  the  AOD  retrieval  process.  Had  the 
ship  Ron  Brown  been  able  to  gather  more  in-situ  upward- 
looking  radiance  measurements  and  collect  some  of  the 
phytoplankton  for  chemical  analysis,  it  might  be  possible 
to  screen  its  reflective  bandwidths  in  the  future.  The 
interrelated  light-reflecting  properties  of  aerosol  and 
ocean  sediment  continue  to  plague  researchers  and  the 
warfighter  interested  in  maritime  safety  during  mine- 
clearing  operations.  For  further  discussion  on  turbidity 
retrievals  and  its  effects  on  remote  sensing  during  the 
same  ACE-Asia  experiment,  see  Rocha  (2001) . 

The  results  of  this  study  of  NPS  AOD  model 
modification  merely  builds  upon  the  foundation  lain  by 

Brown,  1997  and  Smith,  1998.  Further  steps  will  surely 
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incorporate  improved  sensor  technology  such  as  MODIS  aboard 
the  TERRA  satellite  or  the  future  NPOESS. 

Eurther  areas  to  explore  in  regards  to  this  study 
include : 


•  Improved  model  integration  for  vertical  lifting 

and  surface  type  for  dust  source  region,  relative 
humidity  and  subsidence  for  aerosol  growth  and 
water  deposition  (radius  change)  and 

precipitation  from  the  air  column  over  time. 

•  Improved  cloud  mask  ability  in  the  lack  of 
infrared  split -channel  techniques. 

•  Improved  plankton  bloom  detection  and  testing  the 
purely  single-scatter,  specular  reflection, 
spherical  particle  assumptions. 

•  Ground-truth  data  collection  was  hampered  in  this 
experiment  due  to  political  constraints.  Had  the 
ship  Ron  Brown  or  the  C-130  and  King  Air  been 
able  to  spend  more  time  in  the  Yellow  Sea,  a  more 
complete  knowledge  of  the  aerosol  constituents 
would  be  known.  Definitive  data  on  the  size  and 
chemical  makeup  of  this  dust  leaving  the  Chinese 
coast  would  have  allowed  for  further  model 
validation.  The  parameterization  of  the  phase 
function  remains  a  largest  error  source  and  the 
greatest  uncertainty  for  AOD  algorithms.  In  the 
lack  of  in-situ  measurements,  future  improvements 
to  tools  such  as  the  NASA/Goddard  Space  Elight 
Center  MODIS  dust  classification  look  promising. 


54 


•  Add  data  sets  to  this  experimental  series  by 
measuring  AOD  along  with  plankton  bloom-prone 
regions  navigable  for  US  research  vessels, 
including  the  southern  California  coast  and  the 
Gulf  of  Mexico  as  substitutes  for  the  Yellow  Sea 
and  Bohai  Bay.  Currently,  good  results  have  been 
achieved  in  the  vicinity  of  Taiwan;  further  study 
into  the  South  China  Sea  and  the  South  Pacific  is 
warranted . 
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